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The Lyapunov exponents for Anderson localization are studied in a one dimensional disordered 
system. A random Gaussian potential with the power law decay ~ of the correlation function 

is considered. The exponential growth of the moments of the eigenfunctions and their derivative is 
obtained. Positive Lyapunov exponents, which determine the asymptotic growth rate are found. 
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In this paper we consider Anderson localization [l|, y| 
in a one dimensional disordered system with a long-range 
memory. The recent realization of disordered systems by 
using ultra cold atoms [s', '3| in optical lattices and mi- 
crowave realization of the Hofstadter butterfly [^] show 
that the random potential in the experiments are highly 
correlated. The increased interest in the problem of 
Anderson localization in random potentials with long- 
range correlations is also relevant to studies of the metal- 
insulator transition [1, 0] ■ 

Anderson localization in a one dimensional disordered 
system is described in the framework of the eigenvalue 
problem 

e^{x)^-^<P{x)-V{x)<j,{x), (1) 

with a Gaussian random potential V{x). The long-range 
memory of the disorder is modelled by the two point 
correlation function C{x) with the power law decay at 
the large scale 

{V{x')V{x)) = C,{x - x') = , (2) 

\x — x'yi 

where g > 0. It has been shown by various techniques 
that all eigenfunctions are localized for correlated poten- 
tials with the correlation decay rate < g < 1 0, 
Spectral properties of the random operator of Eq. (TO) 
(and its discrete counterpart) were studied [1, [ifl, [Hi- 
Due to the physical interpretation, see discussion in Ref. 
Q , one of the main results is the absence of the absolutely 
continuous spectra for the random Schrodinger operator 
^ with the correlation properties due to Eq. This 
means that the eigenfunctions 4){x) are localized, and 
investigation of Lyapunov exponents is a serious task re- 
lated to localization of the eigenfunctions. 

The Lyapunov exponents are important in spectral 
theory, since they govern the asymptotic behavior of the 
wave functions. They are defined on the asymptotic be- 
havior of the averaged envelope 7s (e) = lima;_^oo ^''^tJ^^^ ■ 
It was shown by rigorous analysis that the positive Lya- 
punov exponents are absent for the absolutely continuous 
spectrum, while the positiveness of the Lyapunov expo- 
nents ensures that the spectrum is pure point [LI, 12]. 

In this paper, we calculate {(j)'^{x)) of solutions of Eq. 
(fT|) for a certain energy e, with given boundary conditions 



at some point, for example (j>{x = 0) and 4>'{x = 0), where 
prime means the derivative with respect to x. Since the 
distribution of random potentials is translationally in- 
variant, it is independent of the choice of the initial point 
as a; = 0. It will be shown that this quantity grows expo- 
nentially with the rate 7(e) ~ lim2:^oo ^"^'^^^^^^ > 0. Note 
that it is different from -7,, , which supposes a knowledge 
of all the even moments [ij, [3, [l^, [la| . 

We develop a general procedure which is suitable for 
calculation of all moments of the wave function and its 
first derivative. To this end the Schrodinger equation ([T]) 
is considered as the Langevin equation and the x coor- 
dinate as a formal time. For the S correlated process 
it can be easily mapped on the Fokker-Planck (diffu- 
sion) equation for the probability distribution function 
V{4>,4>') [H, Unlike this, the two point correlation 
function which corresponds to the stationary process, 
leads to additional integration over the formal "time" 
with a memory kernel. The method of consideration en- 
ables one to observe the exponential growth of ((/)^(.x)) 
with the Lyapunov exponent 7(e) > 0. 

Since the Schrodinger equation ([1]) is a linear stochastic 
equation, equations for the 2n moments of the type 

Mk,i{x) = mx)f[4>'{x)\'), k+l^2n, fc,Z = 0,1,2,... , 

(3) 

can be obtained in the closed form. To this end we rewrite 
Eq. (jlj in the form of the Langevin equation. The x 
coordinate is considered as a formal time on the half axis 
X = T, T € [0, 00) and the new dynamical variables 
m(t) — (f){x), v{t) = u — (j)' {x) are defined. In the new 
variables the Langevin equation reads 

ii^v, V = - \e + V{T)]u, (4) 

where V{t) is now the long-range correlated noise 

It is convenient to set 5 = 1-1-0; and Cq = Ca- In the 
new variables the expectation values of Eq. ^ are now 
Mk^iir) = (m'^z;'). Solutions of Eq. (0]) are obtained as 
functionals 

t;(i) = -/" [e + V{T)]u{T)dT , u{t) ^ [ v{T)dT . (6) 
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Following [l5| we obtain a temporal equation for the mo- 
ments from the Langevin equation Q and its solutions 
(l6|). Differentiating Mk.iij) with respect to r, we obtain 

Mk,i = kMu-u+i - leMk+iA-i - l{V{t)u''+\'-^) . (7) 

The application of the Furutsu-Novikov formula 18] to 
the last term in Eq. ([7]) yields 



{V{t)T[Vm ^lJ^dr'{V{t)V{r')){'-g^ 



(8) 



Here the solution of Eq. ([6]) is used to obtain the func- 
tional derivative of the functional JF[y(T)] — 
Substituting the solution of Eq. ^ in Eq. Q , we obtain 
that the temporal behavior of the moments is described 
by the fractional-differential equation 

Mk,i = kMk^i^i+i - leMk+ij-i + h{l - l)DfMk+2,i-2 , 

(9) 

where the convolution integral in Eq. ([8]) is the fractional 
derivative Df f{t) 



f{r)dT 

{t-Ty+' 



(10) 



Here the correlation function Ca{t) defines the memory 
kernel, or the causal function. Eqs. ([9]) and (fTO|) are rele- 
vant to the fractional Fokker-Planck equations which de- 
scribe a variety of phy sical processes related to fractional 
diffusion [l^, [23, [2l| . An important technique for the 
treatment of the fractional equation is the Laplace trans- 
form. It is worth stressing that both analytical properties 
of this fractional integration and the Laplace transform 
depend on a. 

For — 1 < a < Eq. Q is readily solved by 
means of the Laplace transform. Defining C[Mkj{t)] — 
Mk,i{s), one obtains from Eq. HU]) C[DfMk,i{t)] = 
CaT{—a)s°'Mkd{s), where r(a) is the gamma func- 
tion. For simplicity, disregarding the sign of the cor- 
relation function ([5]), we set Cq, = 2/r(— a). Then, 
we introduce 2n + 1-dimensional vectors M„(t) = 

(^M2n,o, M2„_i4, . . . , Afi,2„-i, Mo,2n) in the "time" 

space and ]VI„(s) — £[M„(t)] in the Laplace space, cor- 
respondingly. Then the solution of Eq. ^ is the Laplace 
inversion of the following vector 



M„(s) 



Anis) 



M„(0) , 



(11) 



where (2n -I- 1) x (2n -I- 1) matrix A„(s) consists of coef- 
ficients from the matrix equation ([9|) . In the limit s —> 
the disorder term of order of s" oo is dominant, and 
the maximal eigenvalues of An can be evaluated at the 
energy e « 0. Following Ref. [1J|, it can be proven that 
for e = the maximal eigenvalues of A„ behaves for large 
n as f2(s) « s"/^(2ri)''/^. Expanding the initial condition 



M„(0) over the eigenfunctions of An, we obtain that the 
maximal growth of the nth moment is 



Mn{t) = t 



-a/3 



-"/3 - (2n)4/3 



Mo(0) . (12) 



The inverse Laplace transform is the definition of the 
Mittag-Leffler function [22]: E^^ai-i{^[2nYlH^-°'l'^^ . 
Asymptotic behavior of the Mittag-Leffler function for 
t — > oo is determined by the exponential function 
exp (2n)'*/(3-a)^ ^ Therefore the exponential growth of 
the nth moment is due to the Lyapunov exponent 

7(0) - (2n)4/(3-") (13) 

for -1 < a < 0. 

For a > the fractional integral diverges. To overcome 
this obstacle, one considers the causal function as a gen- 
eralized function, and a suitable regularization procedure 
can be carried out see e.g., 0, [2lj. Let — 1 < a < iV, 
where A^ > 1 is an integer. Again, using the composi- 
tion rule, one obtains the Riemann-Liouville fractional 
integral ()10|) in the regularized form 



fit) 



DRLfit) 

* f{r)dT 



T{N-a)dt^ Jo (t-r)i+"-^ 
Thus Eq. dSI) reads 



(14) 



Mk,i = kMk+i.i-i - leMk+i.i-i + 1(1 - l)D%LMk+2,i-2 ■ 

(15) 

This fractional equation of the order of a must be 
equipped with A^ — 1 quasi initial conditions: in addi- 
tion to the initial conditions Mfe^;(0), one has to know 
A^ — 1 fractional derivative of Mk,i{T) at r = 0. Applica- 
tion of the Laplace transform to the fractional derivative 
yields [H 

Af-l 

(16) 

In the asymptotic limit s — > we obtain that the solu- 
tion of Eq. (jlSp is approximated by the inverse Laplace 
transform of the vector 



M„(s) 



1 



S - An{s) 



[J„ - BnD"^l^]Mn{0) , (17) 



where /„ is an unit matrix and matrix _B„ consists of the 
off diagonal elements which produce Mk+2 1-2 terms in 
Eq. (HSl). 

Since we are seeking the maximal growth rate of the 
solution of Eq. (jlSp and the initial condition are not 
important for this growth, we choose the initial condi- 
tion as the eigenvector of the maximal eigenvalue of the 
matrix A„. In what follows we consider a temporal be- 
havior of the second moments, described by 3 x 3 matrix 
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Ai (s) . The eigenvalues of the matrix are roots of a cu- 
bic equation [23|. The growth rate is determined by the 
eigenvalue with the largest real part that will be denoted 
by f2. Taking the initial condition in Eq. P71 as the 
eigenfunction of il, namely — Mo(0), wc obtain that the 
dynamics of the second moments is due to the Laplace 
inversion 



see e.g., 0, HH, (H . Therefore, Eq. ^ now reads 



Mi(t) oc£-^ {n{s)-s) 



(18) 



For small s the eigenvalues il(s) correspond to a "weak" 
disorder in the Laplace space. Therefore the high en- 

fy limit is valid w s^/e, where e » s" [H, 
. Substituting this eigenvalue in Eq. ([18]) and ex- 
panding the denominator we have for the integrand 
^;^^oe"-i(l/s)("-i)"+". Carrying out the Laplace in- 
version, we obtain the solution in the form of another 
definition of the Mittag-Leffler function (see e.g., 0jl23|) 

i;„_i.„(et"-i) = Er=o ■ Therefore 



r(nQ: — n+o) 



Mi{t) oc er-i£;„_i,„(er' 



(19) 



Since the argument of the Mittag-Leffler function is pos- 
itive et°'~^ > 0, then the asymptotic behavior is approx- 
imately Eg sjz) ~ z''-^~^''/^ exp{z^^^) for z — > (X) for all 
values 5 ^^,UM- Therefore, when et"^^ — + oo the expo- 
nential growth of the second moment 

Mi(i) (X exp[7(e)i] (20) 

is approximated by the Lyapunov exponent 

7(e) - ei/("-i) . (21) 

Another way to obtain the Lyapunov exponents avoid- 
ing the difficulties related to the iV — 1 quasi initial con- 
ditions in Eqs. ([T5|) and (|16p is to discard the causality 
principle and extend the consideration of the random pro- 
cess on the entire x axis x e (— oo, -l-oo). For this formal 
consideration, the Furutsu-Novikov formula in Eq. ([8]) 
reads 



{I - 1)C„ 



M, 



'J-k+2J~2 



iv) 



[x - yY+^ 



dy for X > , 



f^^rf^ for.<0. (22) 

Setting again Cq — r(— a), we obtain that Eq. (|22p is 
the definition of the Riesz/Weyl fractional derivative 



— Mfe,i = kMk-i.i+i' 
ax 



-leMk 



(23) 

A specific property that we use is the fractional differenti- 
ation of an exponential W" exp(7a;) = 7" exp(7a;). Sub- 
stituting this in Eq. ([7]), one seeks the solution for the 
maximal moment growth Mk^i{x) = exp(±7a;)Mfc^; (a; = 
0), where plus stays for x > and minus for a; < 0, re- 
spectively. One readily checks that the both cases yield 
the same algebraic equation 



7M„ = An{-/)Mn, 



(24) 



where the moment vector M„ is defined above and the 
matrix Anij) is defined from Eq. Therefore, 
$7(7) = 7"/e, where conditions 7 <C e and 7" ^ e are 
used. Solutions of Eq. (|23)) for 7(e) coincide exactly with 
the ones obtained in Eqs. (fT^)) and (|^T|) for all values of 



This solution for 7 also yields conditions of validity of 
the solution (|2T|) for different values of energy e. Indeed, 
for < a < 1 Eqs. ^ and ([M]) describe an 

exponential growth for asymptotically large energies e ^ 
1, since, in this case, 7(e) ^ e when e 3> 1. On the 
contrary, when a > 1 the solution of Eq. ([2T|) is valid 
for e ^ 1. This follows from the condition 7" <C e. Note 
that for large negative values of the energy Q 2A/[e[, 
what corresponds to a simple pole in Eq. (|18p . and this 
is just the Lyapunov exponent 7(e) ~ 2-^/[ef . 

In conclusion, we studied the Lyapunov exponents for 
Anderson localization in a one-dimensional disordered 
system with a long-range memory. The averaged behav- 
ior of the second moment of the eigenfunction is calcu- 
lated, and its asymptotic exponential growth for \x\ — > 00 
is determined by the Lyapunov exponents for different 
values of the energy e. The main result of the study is the 
existence of the positive Lyapunov exponents 7(e) > for 
the rate q = l + a>Ooi the power law decay of the cor- 
relation function. It is relevant to the exponential local- 
ization of the eigcnfunctions of the random Schrodinger 
operator of Eq. ([1]). 
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